arXiv:1506.05656vl [astro-ph.CO] 18 Jun 2015 


Helium reionization in the presence of self-annihilating clumpy dark matter 

Bidisha BandyopadhyajQ 

Department of Physics and Astrophysics, University of Delhi, Delhi, India 
Dominik R.G. SchleicheiQ 

Departamento de Astronomia, Facultad Ciencias Fi'sicas y Matemncas, Universidad de Concepcion, 

Av. Esteban hurra s/n Barrio Universitario, Casilla 160-C, Concepcion, Chile 

The reionization of helium describes the transition from its singly ionized state to a doubly-ionized state in 
the intergalactic medium (IGM). This process is important for the thermal evolution of the IGM and influences 
the mean free path of photons with energies above 54.4 eV. While it is well-known that helium reionization is 
mostly driven by the contribution of energetic quasars at z < 6, we study here how helium reionization proceeds 
if there is an additional contribution due to the annihilation of dark matter. We explore the effects of different 
dark matter profiles for the dark matter clumping factor, which can significantly enhance the annihilation rate 
at late times. We find that the presence of dark matter annihilation enhances the He ++ abundance at early 
stages where it would be zero within the standard model, and it can further increase during structure formation, 
reflecting the increase of the dark matter clumping factor. The latter is, however, degenerate with the build-up 
of the quasar contribution, and we therefore expect no significant changes at late times. We expect that future 
studies of the He + Lyman a forest may help to assess whether the evolution is consistent with the contribution 
from quasars alone, or if an additional component may be required. 


I. INTRODUCTION 


About 24% of the baryonic mass in the Universe con¬ 
sists of primordial helium, implying that the ionization 
state of helium plays an important role in the thermal 
evolution of the intergalactic medium (IGM) which in¬ 
fluences the mean free path of photons with energies 
greater than 54.4 eV. Measurements of the He + Lya op¬ 
tical depth suggest that helium has been reionized and 
now occurs predominantly in the double-ionized state. 
These measurements are based on Hubble Space Tele¬ 
scope (HST) observations of Q0302-003 |f39l l44h . ob¬ 
servations from the Hopkins Ultraviolet Telescope and 
the Far Ultraviolet Spectroscopic Explorer (FUSE) of HS 
1700+6416 [[Mfl . HST and FUSE observations of HE 
23474342 MBolMlbl, and HST observations of PKS 
1935-692 S7 sdSS J2346-0016 JIMl, Q1157+3143 
HI and SDSS J1711+6052 0], spanning a redshift 
range 2.2 < z < 3.8. 

The reionization of helium (He + —> He ++ ) requires 
highly energetic photons with energies above at least 
54.4 eV, and is therefore accompanied by a significant 
increase in the thermal energy of the IGM 11341 [38 1. 
Such changes in the thermal evolution can have a 
significant impact on structure formation, in particular 
in the context of the formation of lower-mass objects 
like dwarf galaxies lf55ll . Understanding the evolution of 
helium reionization is therefore important to understand 
the evolution of the IGM as well as the formation of 
galaxies within the IGM. 


Observations from galaxy surveys and CMB signals 
suggest that most of the matter content of the universe 
is non-baryonic in nature. In particular, there is a pres¬ 
sureless component representing about 31% of our Uni¬ 
verse which is called dark matter (from the PLANCK 


* bidisha@physics.du.ac.in 
t dschleicher@astro-udec.cl 


survey |5il), which leaves an imprint on various astro- 
physical and cosmological observations only due to its 
gravitational interaction. Many particle physics models 
suggest that the dark matter consists of weakly interact¬ 
ing massive particles (WIMPs) with typical masses of 
10 — 1000 GeV which may annihilate into photons, elec¬ 
trons or positrons and neutrinos. 

We assume here an approximately equal distribution 
among the potential annihilation products, as expected 
from generic dark matter models. We however note that 
an annihilation signal would not be visible if the annihi¬ 
lation products consist of neutrinos only. As a result, an 
excess of cosmic ray signals from the center of galaxies 
is expected due to the self-annihilation of these dark 
matter particles. From recent studies dark matter masses 
M x < 26 GeV are ruled out for a thermally averaged 
annihilation cross-section (<r x v) = 3.0 x 10” 26 cm 3 s _1 
and perfect absorbtion or for a more realistic absorbtion 
M x < 5 GeV ll49l l. The dark matter candidates which 
are massive and have a negligible free-streaming length 
compared to the first protogalaxies are considered to 
be ’cold’ while the warm dark matter particles (WDM) 
are comparably less massive and their free-streaming 
length is comparable to the size of protogalaxies. WDM 
is a possible alternative to CDM since it can explain 
the presence of cuspy halos in our universe Emnim 
Light neutralinos of some SUSY models can be CDM 
candidates, while sterile neutrinos j|j, [5, HU], majorons 
IH M SI] an d light dark matter particles (LDM) fl2h 
are proposed to be among the candidates for WDM. 
While WDM can explain cuspy halos in our Universe, it 
is ruled out by recent CMB observations and other data 
sets B We refer in particular to recent reviews §0 
for the different dark matter models. 

If dark matter consists of self-annihilating particles, 
one may expect to observe potential annihilation signals 
in regions where the dark matter density is particularly 
high, as in the center of our own galaxy. For instance, 
strong X-ray emission in the 511 keV line from our 
galactic center can be interpreted as a sign of self- 
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annihilating dark matter 11131. l45l|70l[ . Other observations 
suggest an excess of GeV photons [ 14], microwave pho¬ 
tons 0 and positrons (23:] in the galactic center. An 
interpretation in terms of self-annihilating dark matter 
is favored by the approximately isotropic emission of 
such particles and the absence of correlations with the 
Galactic disk. 

Apart from the search for dark matter signals in 
our galactic center, models for their impact on various 
cosmological and astrophysical processes have been 
suggested in the literature. Both particle physics and 
cosmological data provide bounds on their annihilation 
cross section SSIMSSSEilSIil and 
potential signatures in the recombination spectrum have 
been explored by Chluba (2009). The annihilation of 
dark matter particles leads to the production of energy, 
part of which is injected into the IGM, where it can 
affect the ionization and thermal history of the Universe. 
In general the ionization of hydrogen and helium 
occurs through similar processes, and can therefore be 
described with a similar framework (e.g. 1,35!]). Since 
dark matter particles are massive, their annihilation 
produces photons which have very high energies and 
typically provide minor contributions to the ionization 
of hydrogen and helium, unless extreme scenarios 
are considered msa. On the other hand more 
energetic photons are required to ionize He + to He ++ 
and such photons are still rare among the first sources of 
light. Hence the contribution from the self-annihilation 
of dark matter is potentially more relevant in this case 
and will be explored in greater detail in this paper. 

This paper is structured as follows. In 3U] we out¬ 
line the different sources which affect the He ++ frac¬ 
tion. In subsection 311 A1 we explain the ionization equa¬ 
tion, and in subsection 311 B1 we describe the contribu¬ 
tion of quasars to the formation of He ++ . In 311 Cl we 
discuss the effects of dark matter annhilation assuming 
a uniform dark matter density, while in 311 Dl we take 
into account the clumped dark matter density taking into 
account different halo profiles. In sections 31111 we de¬ 
scribe our main results, and discuss our conclusions in 
3IVI The annhilation rates are evaluated in detail in the 
APPENDIX 


II. HELIUM REIONIZATION WITH QUASARS AND 
DARK MATTER 

Singly ionized helium with a ionization threshold of 
-Ejte+.ion = 24.6 eV is produced along with the ioniza¬ 
tion of atomic hydrogen (E’n+.ion = 13.6 eV), which 
is generally considered to occur as a result of photoion¬ 
ization due to radiative feedback from the first stars and 
galaxies 0 22,, 64-J. He + has a similar recombination 
rate as H + , while the recombination rate is enhanced for 
He ++ by a factor of 5. In studies of helium reionization, 
it is thus typically assumed that the single-ionization of 
helium is more or less co-eval with hydrogen reioniza¬ 
tion, while the ionization of He 1 + occurs at a later stage. 
The ionization of He + to He ++ further requires photons 
of higher energy owing to a greater ionization thresh¬ 


old (-EHe++,ion = 54.4eU). Even in models consider¬ 
ing reionization due to massive Pop. Ill stars, the num¬ 
ber of ionizing photons at 54.4 eV is suppressed by at 
least a factor of 20 compared to the peak of the emis¬ 
sion. However, observational studies suggest that reion¬ 
ization is predominantly driven by low-luminosity galax¬ 
ies at high redshift |]lg], implying a more conventional 
stellar population. As a result, the expected drop in the 
number of He + -ionizing photons is even more signifi¬ 
cant, requiring the emission in quasars or potentially due 
to the self-annihilation of dark matter. 

In the following subsections, we outline the frame¬ 
work to describe how helium is doubly-ionized from its 
singly ionized state due to photons produced in quasars 
or as a result of dark matter annihilation. As shown in 
our calculation below, the formation of He ++ occurs dur¬ 
ing epochs where helium is already singly-ionized. We 
therefore adopt the single-ionization of helium as a gen¬ 
eral assumption in our calculation. While this could in 
principle lead to an overestimate of the He ++ produc¬ 
tion in epochs where the helium is still neutral, our re¬ 
sults show that the latter happens to be inefficient dur¬ 
ing such times, therefore providing an a posteriori jus¬ 
tification. For this study, we adopt a generic set of 
cosmological parameters approximately consistent with 
the recent PLANCK lf56ll survey with fl m = 0.31 , 
fl A = 0.68, fl b = 0.033, H = 100 h km s^Mpc” 1 , 
h = 0.6711, n = 0.96, as = 0.8 and the critical density 
Pcrit = 3.1 x lO 17 /!" 1 M 0 Mpc -3 . 


A. IONIZATION EQUATIONS 

While the main focus of this investigation concerns 
the epoch of helium reionization, a residual degree of 
doubly-ionized helium can be produced already in the 
epoch of recombination. In the absence of dark matter 
annihilation, it is well-known that helium recombines 
completely and becomes fully neutral. However, in the 
presence of dark matter annihilation, there will be an 
additional ionization term which prevents the He ++ 
abundances from recombining to zero, and instead a 
new equilibrium will be established corresponding to a 
non-zero doubly-ionized ionization degree. As already 
shown by Chluba (2009), helium is in the singly-ionized 
state until z ~ 2000. Subsequently, a residual ionization 
of ~ 10 -4 is maintained in the singly-ionized state. 
While the recombination rates of He ++ are higher by 
about a factor of 5 compared to He+, residual abun¬ 
dances of up to ~ 10 -5 can nevertheless be maintained. 

As a minimal model to obtain the residual ionization 
degree from the recombination epoch, we will in the fol¬ 
lowing consider the ionization due to dark matter anni¬ 
hilation, the recombination to and the collisional ioniza¬ 
tion of He + . In addition to the processes outlined above, 
we consider the double-ionization of helium by energetic 
photons from high-redshift quasars. As we will show be¬ 
low, these are likely to give the dominant contribution 
leading to complete helium reionization in the intergalac- 
tic medium. The rate at which the mean ionization frac- 
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tion of He±± evolves is then given by 


*^He++ 

dt 


kQSO ± 3“ (1- *£_fTe++ )t^e/5//e+ 

C'a J 4(T)n e 5_; ?e ++. (1) 


Here, the first term describes the ionization rate due to 
quasar photons, the second term describes the contribu¬ 
tion of dark matter, the third term describes the colli- 
sional ionization by thermal electrons and the fourth term 
describes radiative recombination. The fraction of He ++ 
is expressed as i He ++ = ”^ ++ , while the number den¬ 
sity of thermal electrons is n e , and /3fj e + id is the colli- 
sional ionization coefficient, which evolves with the mat¬ 
ter temperature and is given as 


ionization rate of helium due to ionizing photons from 
quasars is given as 


kQSO 



N dct> 
n,H e dLs 


N = 2 x 10 55 


Lb 


10 12 L 


o 


(5) 

( 6 ) 


where N is the total number of photons emitted per sec¬ 
ond by a quasar which have energy greater than the ion¬ 
ization threshold of He ++ . Here dd/dLn is the quasar 
luminosity function (QLF), which is defined as the num¬ 
ber of quasars per unit volume per unit luminosity in¬ 
terval. Following the work of Hopkins et. al. lidll . we 
obtain the following expression: 


P He+ = 5.68 X l 0 -12 T l/2 e -631515/T 

cm 3 s -1 . ( 2 ) 

The recombination coefficient a A (T) d is given as 
a A (T) = 3.36 x io- 10 T " 1/2 



1 + 


10 5 I 


1/2 


X 



cm 3 s 1 


(3) 


and C = (p 2 ) / (p) 2 is the baryonic clumping factor 
which we describe as 

f 1 (z > 15) 

C = l 1 + (6 < z < 15) (4) 

(3 (z < 6 ). 

In this context, we note in particular that a uniform 
baryon density is assumed at z > 15, before a significant 
amount of structure formation has occured, and a generic 
value of C = 3 is adopted at 2 < 6 , as also employed in 
other studies (e.g.[34j|). A linear interpolation is pursued 
in between. We note that the baryon clumping factor 
used here is smaller than in studies of hydrogen reion¬ 
ization EH, as helium reionization is strongly driven 
by rare objects like quasars and therefore proceeds on 
large spatial scales. 


The processes outlined here are important at different 
cosmological epochs, and will lead to characteristic sig¬ 
natures of the evolution of He ++ at different cosmolog¬ 
ical times. In particular the evolution of the recombina¬ 
tion and collisional ionization rate depends on the tem¬ 
perature of the gas. We have used the RECFAST code 
1631] to follow the gas temperature up to z = 15 and be¬ 
low z = 15 we assume the IGM temperature rises to 
T = 10000 K during cosmic reionization. 


B. CONTRIBUTION FROM QUASARS 

Quasars and AGNs produce hard photons, which sub¬ 
stantially contribute to the production of doubly-ionized 
helium. Following Furlanetto et. al. 0 , the doubly 


dtp _ _ <fi* _ 

dlogL B ~ (£/T *) 71 +(L/TQt' 2 ’ 

where </>* is the normalization, L* is the break luminos¬ 
ity, 71 is the faint-end slope and 72 is the bright-end 
slope. 

The equation for the QLF corresponds to a double 
power law which is also evident from Fig. iflall . The pa¬ 
rameters involved in the expression for the QLF evolve 
as 


log L* 

71 

72 


€ 


(log T*)o + + ^L3? 3 , 

( 7 i)o x 10 fe ^«, 

_ 2 ( 72)0 _ 

2Q fc 72,l? _|_ 10 fc 72,2f ’ 


log 


f 1 + z 

\i + 2/ 


( 8 ) 


where kj, 1 , k[/ 2 and k ^3 are free parameters and Zf = 2. 
The best fitparameters for the above model are given in 

[Table □ 0 ], 


TABLE I: Best fit Parameters for QLFs lf4lTl 


Parameters 

Values 

logtj> * 

—4.825 ± 0.060 

logLt 

13.036 ± 0.043 

kLi 

0.632 ± 0.077 

kL2 

-11.76 ±0.38 

kL3 

-14.25 ±0.80 

(71)0 

0.417 ±0.055 

k- yi 

-0.623 ±0.132 

(72)0 

2.174 ±0.055 

k l2 i 

1.460 ± 0.096 

^ 72,2 

-0.793 ± 0.057 


where </>* is in units of Mpc -3 and L, is expressed 
in terms of Lq. The effect of ionization by quasars 
has been included from z ~ 6 since that is the epoch 
where the first quasars have been observed, and from 
which we expect the adopted quasar luminosity function 
to be valid. The ionizing photons from quasars lead to 
a rapid double-ionization of helium. The timescale for 
this process is much smaller than the timescale for colli¬ 
sional ionization. Figure |[Tal shows the QLF as a func¬ 
tion of luminosity for different redshifts. The impact of 
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(a) 


(b) 


FIG. 1: (a) The quasar luminosity function dLog [ / / LQ at different redshifts. 

(b) Evolution of the He ++ fraction in the presence of quasars without dark matter annihilation. 


quasars on helium reionization is given in Figure ED. 
which shows the evolution of the He ++ fraction as a 
function of redshift assuming different clumping factors 
C for the baryons. Due to the effect of recombination, 
the amount of He ++ decreases with the baryonic clump¬ 
ing factor. While the calculation here implicitly assumes 
that all quasar spectra are the same, observational investi¬ 
gations have shown that characteristic quasar spectra can 
be obtained when considering a larger sample of sources, 
so that a well-defined statistical average exists ll 6 lU 62 | . 
We expect that the latter will be sufficient to describe the 
large-scale average of helium reionization. So far, even 
the highest-redshift quasar known to date at z=7.085 still 
shows the same characteristic features as quasars at lower 
redshift |52i], therefore providing no evidence for a cos¬ 
mological evolution. If such an evolution is found in the 
future, it could further modify the results for the cosmic 
reionization of helium obtained here. 


C. EFFECTS OF DARK MATTER 

Dark matter constitutes about 28% of our Universe 
[5bi|• While dark matter candidates need to be extremely 
stable with long lifetimes and small annihilation cross 
sections, the annihilation of even a small fraction of the 
total abundance can have a significant effect on the reion¬ 
ization of hydrogen and helium. We adopt here the typ¬ 
ical framework employed to consider the impact of dark 
matter on hydrogen reionization |64j], and adapt it for he¬ 
lium. In case of a uniform dark matter distribution, we 
obtain 

k D M = mXi(<rv) ( E m ) 

\ HeTpon / 

n 2 sf i (fe )( 1 + ^) 3 > ( 9 ) 


where (av) = 3 x 10 ~ 26 cm 3 s~ 1 is the thermally av¬ 
eraged cross section for annihilation, Xi is the fraction 
of energy that goes into ionizing He + and is defined as 
Xi = E/E p h where E is the energy that goes into ion¬ 
ization while E p h is the energy of the incoming particle. 
The total energy required for ionization is E = N x E lnn 
where N is the total number of ionization and E ion is the 
ionization energy. From [27] we have the mean energy 
per ion pair W (i.e. W = E p h/N ) which is parametri¬ 
cally given as W = Wo (1 + Cx° ). Thus substituting, we 
have Xi = Ei on /W. For photons with energies greater 
than 1 keV and a cosmological mixture of hydrogen and 
helium gases the value of the constant parameters are 
W 0 = 16400 eV, C = 11.7 and a = 1.05. For He+ 
ionization we assume 772 = 8/0.24, while the ionization 
energy £ H e++, ion = 54.4 eV. 

Since the dark matter annihilation rate depends on the 
dark matter density, we start our calculation at z = 1000 
to account for the formation of the residual amount of 
He ++ from the annihilation at early stages. We will in 
the following consider both warm dark matter particles 
with masses of about 20 MeV, and cold dark matter par¬ 
ticles with masses of about 80 GeV. 


D. DARK MATTER CLUMPING 

In the previous section we have assumed that the 
dark matter is uniformly distributed. While this uniform 
density of dark matter can be a good approximation 
at early stages, it is expected to clump considerably 
during structure formation. This non-uniform density 
distribution can be accounted for with a clumping factor 
for the dark matter, which effectively describes how the 
dark matter annihilation is enhanced due to the dark 


x 
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matter clumping or by methods as mentioned in j43f |. 
The latter depends both on the abundance of dark matter 
halos, as well as their characteristic density profiles. 

Here, we consider the Moore, NFW and Burkert pro¬ 
files for the dark matter [ 18, 51, 53l |54l 0 and adopt 
the approach by Cumberbatch et al. |26j] to calculate the 
clumping factor. We first consider the total annihilation 
rate resulting from dark matter halos with a minimum 
mass M m i n and a maximum mass M max at a redshift z 
which is given as 


r halos(z) = (1 + z) 3 

rMmax 1 

X J M dM—(M,z)R(M,z),m 

where dn/dM is the mass function which gives the num¬ 
ber of halos with mass M and R(M , z) is the rate of 
annihilation inside a halo of mass M at the redshift z. 
As discussed in the appendix, we will typically adopt a 
lower limit of 10 6 M 0 for the halo masses, even though 
it can be as low as 10 -12 M 0 in WDM models [|26[|. As 
shown in Fig. [5] even significant changes in this parame¬ 
ter do not strongly affect our results. 

We evaluate the halo mass function using the Press- 
Schechter formalism ill. The annihilation rate inside a 
halo of mass M at redshift z is then 

R{M,z) = P 2 {r)^r 2 dr (11) 

2 J r= o 

where r V i r is the virial radius. The annihilation rate in 
the smooth background is 


Fsmooth{z) = \ ^2 ^ P X ( Z ) (12) 

z m x 

and thus the clumping factor at a given redshift(z) is 
given as 


C halo (z) = l+ ^ halos{z) (13) 

-t smooth) 

The annihilation rate inside the halos can now further 
be determined for different dark matter density profiles. 
In general, the density distribution can be described as a 
continous function of r (26']: 

P ^ ~ (r/r s ) 7 [l + (r/r s )“](/ ? -T')/“' (14) 

We will focus in particular on the NFW, Moore and 
Burkert profile 10 BTl I53L El 0 . The derivation of 
the rates of annihilation for each of the profiles is given 
in the Appendix. 

The ionization rate from dark matter annihilation 0 
is now modified to include the clumping factor and can 
be written as 


&DM = mXiChalo{z)(crv) 



III. RESULTS 


The evolution of helium reionization in the presence of 
quasars and in the absence of dark matter annihilation is 
given in Fig. m for baryonic clumping factors ranging 
from 0 (underdense regions with inefficient recombina¬ 
tion) to 3 (overdense regions with enhanced recombina¬ 
tion). The results show a very steep increase of the He ++ 
fraction as a function of redshift, with helium reioniza¬ 
tion being completed at z ~ 2.5. The latter reflects both 
the shift of the quasar luminosity function towards higher 
luminosities, as well as the increase of the quasar abun¬ 
dances at a given luminosity in Fig. HU ■ The results 
depend only weakly on the baryonic clumping factor in 
the calculation. 

The rate of dark matter annihilation depends on the 
number density of the dark matter particles given by 
n x = Pdm/trt-^, where pum — - k/mPcrit 5 with p cr it the 
critical mass density and ttdm ~ 28%. The dark matter 
content of our Universe may consist of massive cold dark 
matter particles ( m x > 10 GeV) or warm dark matter 
with characteristic masses m x ~ MeV. Due to the large 
difference in the masses of warm and cold dark matter, 
there is a corresponding difference in the number density 
of dark matter particles for these two models. In Fig. 0 . 
we compare the results when considering the contribu¬ 
tions of cold and warm dark matter to helium reioniza¬ 
tion, finding that the contribution from warm dark matter 
is clearly enhanced due to the higher annihilation rate. 
However, the warm dark matter scenario is ruled out by 
the PLANCK survey 10 . so we will subsequently focus 
on the cold dark matter scenario. 

We can generally distinguish a number of character¬ 
istic phases in the evolution of He ++ : At z < 1000, 
the He ++ abundance decreases but remains non-zero, re¬ 
flecting the decrease of the annihilation rate and the de¬ 
crease in the gas temperature, which is initially coupled 
to the temperature of the CMB. At z ~ 30, the effects 
of dark matter clumping become noticable which starts 
increasing the dark matter annihilation rate and hence 
the He ++ fraction, which is further enhanced when the 
gas temperature increases to ~ 10 4 K due to hydrogen 
reionization. Subsequently, we find a further moderate 
increase until the contribution from quasars kicks in at 
z = 6. 

We notice that these phases can be identified both in 
the calculations with cold and warm dark matter and 
for all the halo profiles considered, including the NFW, 
Burkert and Moore profile. The evolution of the dark 
matter clumping factor is given for these cases in Fig. 0. 
As expected, we find the highest clumping factor for the 
Moore profile, which is the steepest density profile, while 
it is lowest for the Burkert profile, in agreement with the 
study by Cumberbatch et al. lf26fl . We expect that the real 
dark matter clumping factor will lie in between these ex¬ 
treme cases. 

A comparison of the implications of the different 
clumping factors for the cold dark matter scenario in¬ 
cluding also a case with a uniform dark matter distribu¬ 
tion is shown in Fig. 0. It is clearly visible that the 
He ++ abundance is reduced by at least a few orders of 
magnitude in case of a uniform DM density. Apart from 
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FIG. 2: Evolution of the dark matter clumping factor for the NFW, Moore and Burkert dark matter density profiles. 





(a) 


(b) 


(c) 


FIG. 3: Evolution of He ++ in the presence of quasars and dark matter annihilation for a warm and cold dark matter 
scenario and compared with the case when there is no dark matter annihilation. ll3al Burkert profile. Obi NFW 

profile and ED Moore profile. 


the increase in the He ++ fraction due to the temperature 
increase at z ~ 15, the evolution of the ionized fraction 
is generally decreasing with decreasing redshift at 2 > 6 , 
reflecting the decreasing density of dark matter and the 
decreasing annihilation rate. This behavior is inverted 
when clumpy dark matter is considered, as the dumpi¬ 
ness increases with redshift and compensates for the de¬ 
crease of the mean dark matter density. As a result, the 
contribution of dark matter annihilation may rise rapidly 
at z < 6 , when also the abundance of luminous quasars 
is rapidly increasing. At 2 < 6 , the fraction of He ++ 
increases both due to the increased temperature of the 
IGM, as well as due to the increasing dark matter clump¬ 
ing factor, which varies between ~ 100 and ~ 3000 at 
2 = 6 for the scenarios considered here. As a result, the 
He ++ abundance is clearly enhanced for clumpy dark 
matter compared to a uniform distribution, and thus rises 


more steeply. A close zoom-in shows that the rise is even 
steeper for the Moore profile compared to NFW or Burk¬ 
ert, but then converges from 2 ~ 4 when the abundance 
from quasars becomes important. The completion of he¬ 
lium reionization may thus change slightly by A 2 ~ 0.5 
due to the additional contribution, even though the lat¬ 
ter is hard to distinguish from the expected evolution in 
the presence of quasars. As the formation of quasars is 
expected to follow the build-up of structures, we expect 
that both contributions may evolve in a similar fashion, 
thus preventing a significant change in the evolution at 
late times. 

Finally, we also checked how the results depend on the 
adopted integration range when we determine the dark 
matter clumping factor. As visible in Fig. 0, there are 
mostly minor changes by at most one order of magnitude, 
even in cases where the integration range is changed by 
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FIG. 4: Evolution of the He ++ fraction for a cold dark matter scenario in case of a uniform distribution, and for the 
case of clumpy dark matter assuming Moore, NFW and Burkert profiles. The mass range of the dark matter halos is 
considered to range from 10 6 M 0 to 10 12 M 0 . It also compares the ionization fraction when no dark matter 

annihilation takes place. 


several orders of magnitude. Even though there is a phys¬ 
ical uncertainty on the appropriate range to consider, we 
find that the overall behavior is rather robust and we can 
again identify the same characteristic phases in these cal¬ 
culations. 


IV. DISCUSSION AND CONCLUSION 

We have investigated the impact of self-annihilating 
clumpy dark matter on the evolution of the He ++ abun¬ 
dance over cosmic history. We find that the contribution 
of dark matter is potentially relevant before quasar feed¬ 
back becomes important, and provides a residual He 1 ^ 
fraction already at high redshift that depends on the dark 
matter model. 

In models with a uniform dark matter distribution, 
the contribution from dark matter annihilation decreases 
with decreasing redshift due to the decreasing dark mat¬ 
ter densities, which is generally reflected in decreasing 
abundances of He ++ before the feedback from quasars 
becomes more relevant. In models of clumpy dark mat¬ 
ter, the dumpiness however increases with time as a re¬ 
sult of structure formation, and therefore both the anni¬ 
hilation rate and the He ++ abundance increases with de¬ 
creasing redshift. We further note that the thermal evo¬ 
lution of the IGM is relevant for the residual abundance 
of He++, as the increased temperatures due to hydro¬ 
gen reionization give rise to reduced He ++ recombina¬ 
tion rates and therefore a higher net abundance. A simi¬ 
lar increase occurs at z < 6 when the temperature of the 
IGM raises to ~ 20000 K. 

While in the case of a uniform dark matter distribu¬ 
tion, the annihilation rate eventually becomes negligible 
at low redshift, there is a considerable enhancement in 


the presence of clumpy dark matter, corresponding to a 
clumping factor of 10 2 — 3 x 10 3 at z = 6 and a clump¬ 
ing factor of 3 x 10 3 — 3 x 10 5 at z = 1, which out¬ 
weighs the decrease of the dark matter density with red¬ 
shift. As a result, there can be a substantial contribution 
at z ~ 3 where helium reionization is expected to ap¬ 
proach completion. The effect is however highly degen¬ 
erate with the contribution from quasars, as both their 
abundance and peak luminosity is expected to increase 
towards z ~ 2 — 3, reflecting the formation and further 
growth of cosmological structures. We therefore expect 
that such a contribution may be present, even though it 
should be difficult to infer how much the completion of 
helium reionization has changed as a result. 

To determine such potential contributions, it is thus 
necessary to both accurately map out the evolution of he¬ 
lium reionization, extending current observational stud¬ 
ies i, HI HU M B 51 El .59, 66, 68. 5! m, but 
also to pursue a precise number count of X-ray sources 
throughout the evolution of the Universe, including both 
Compton-thin j65|| and Compton-thick |69J sources. 

Rather than considering the completion of helium 
reionization, it may be helpful to probe the cosmic evolu¬ 
tion of He ++ at earlier stages by extending current stud¬ 
ies of the He + Lyman a forest and comparing the latter 
with the characteristic evolution of the quasar population 
over time. It is at least conceivable that the slope of the 
evolution provides additional information, which will be 
valuable to distinguish different scenarios. A further as¬ 
sessment of the uncertainties in the observed population 
and its implications for the available amount of X-rays 
for helium reionization will be necessary in this respect. 

We note here that our calculations adopted a mass of 
80 GeV for the dark matter particles, along with a cross- 
section of (av) = 3 x 10~ 26 cm 3 s _1 as consistent with 















(a) 




(b) (c) 

FIG. 5: Evolution of the He 4 1 fraction for cold dark matter assuming [a] a Burkert profile, [b] a Moore profile and 
[c] an NFW profile. We show that the results depend only weakly on the mass range considered, where 

Mi -A 10 6 - 10 12 M 0 , M 2 -a 10 4 - 10 12 M 0 , M 3 -a 10 " 12 - 10 6 M 0 and M 4 -A 10 " 12 - 10 6 M 0 with 

M cut = 10 6 M 0 


the Planck data {56]. We note that for larger cross sec- created, 
tions, the annihilation rate would be enhanced, and the 
contribution to helium reionization thus more significant, 
while the effect becomes less important for lower cross 
sections. Similarly, the mass of the dark matter particles 
determines their number density. In such a case, we still 
expect the same effects to occur, as the dark matter ef¬ 
fects considered here depend only on the ratio (av)/m x 
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The NFW profile is a special case of eqn. dEfli with the 
values a = 1, /3 = 3 and 7 = 1 . The rate of annihilation 
f?(M, z) inside such a halo can be calculated to give 


R(M, z) 


1 {ayv)_ 2 47 t / r vir (z,M) 

2 Ps 3 \c vir {z,M) 


x 



3 

( 1 . 1 ) 


[1 + c vir {M,z)] 3 
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where 


ten as a continous function of the form 


Ps(M, z ) 



X 


1 


p{r) 


_ Ps _ 

[1 + (r/r s )][l + (r/r s ) 2 ]' 


(3.1) 


For such a profile, the annihilation rate follows as 


2. MOORE PROFILE 


For the Moore profile, we have a = 1.5, j3 = 3 and 
7 = 1.5. The annihilation rate inside such a halo is cal¬ 
culated as 


u( M r \ _ 1 (°x v ) n 2 4n ( r V ir(z,M) V 
nyivl,^) — 2 m 2 Ps 3 ) 

R Tl (C-vir t ^min) 


where 


Tl {C-viri %min) — 3 (1 -j-Xmir.) 2 


+ ~ 


l °9 


1 

1.5 


°(l+a£L) 

,L(l + 4ir) 

1 1 


l + c. 


T3" 


l-\-x 


TS~ 


( 2 . 1 ) 


( 2 . 2 ) 


and 


x 


Ps(M, z ) 




+ 5 


1 


-1 


where x m i n = 10 


3. BURKERT PROFILE 

Apart from the NFW or Moore profiles, one often con¬ 
siders the Burkert density profile, which can also be writ- 


X(M,z) = ±te 


^ 2 2 tt 


Itp : 


( r vir (z,M) Y 

\ c„ ir (z,M) J 


1 2(l+c„i r ) 2 tan 1 ( c vir) 2(1+c^ r ) 


(3.2) 


where 


p s (M, z) = 


M 


47T 


( c vir ) 


T X 


[log{ 1 + (Air) + \log{ 1 + c 2 ir ) - tan 1 11 {c vir )\ (^.3) 


In the above expressions c V i r is the concentration 
parameter and is defined as c V i r (M,z) = K . We 
can evaluate the collapse redshift z c with the condition 
(1 + z c ) = where d sc (z) = cr(M*(z)) with 

M* = 0.015 M. The virial radius rw is defined as 

r.i r - U* , I„ order ,o clcula.e ,he value 
of Ml for very small halos we choose a cutoff mass 
M cut below which the concentration parameter remains 
invariant. As 0 , we adopt a typical mass scale of 
10 6 Mq, corresponding to the mass of the halm where 
the first stars are expected to form fl [13 S3. To 
evaluate the mass integral, the mass limits have been 
chosen to be M m i n = 10 6 Mq and M max = 10 12 Mq. 
We have checked that our results are not sensitive to this 
choice. 
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